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We present the first analytical study of phase transitions in ferroelastic-ferroelectric epitaxial thin 
films on exactly solvable model. The emerging domain structure with domains of equal width 
(which may be exponentially large on a "soft" substrate) always remains stable irrespective of the 
film thickness. The dielectric response of an epitaxial film is smaller than that of a free film, in 
striking contrast with assertions in the literature. 
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Formation of ferroelastic domains in thin films due to 
elastic misfit between the film and the substrate was pre- 
dicted by Roytburd E|. The domain patterns are the 
focus of extensive studies, especially in epitaxial films of 
perovskite ferroelectrics that are improper ferroelastics 
|^| . However, the main results in this field are obtained 
by numerical methods or by making seemingly reason- 
able yet unjustifiable approximations. Additionally, all 
previous results refer to the case of equal elastic mod- 
uli of the film and the substrate. As a result, the physics 
of the epitaxial ferroelectric-ferroelastic films remains ob- 
scure, and this has led to generally accepted yet erroneous 
statements in the literature. 

In this paper we consider an exactly solvable case of a 
ferroelectric-proper ferroelastic thin film described by a 
one-component order parameter, which is either a strain 
tensor component or a polarization component |3| ex- 
hibiting a second order phase transition. The model al- 
lows an analytical treatment at all temperatures, includ- 
ing that in the vicinity of the transition. 

We shall begin with the case of the film in the absence 
of an external electric field. The film is assumed to be 
perpendicular to the z-axis and the order parameter the 
u xy component of the strain tensor. It is assumed to 
be attached to an elastically isotropic substrate with the 
shear modulus /i, Fig. 1. The clastic moduli of the fer- 
roelastic are supposed to be the same as in the substrate 
with the exception of the "soft" modulus corresponding 
to the u xy component of strain. Thus, the Landau ther- 
modynamic potential has the form 

F = J dV[2Au 2 xy + 2D(Vu xy f + Bu xy 

+M («? fc - 2*4)] (1) 

where A = a(T—T c ), with a,D,fj, positive constants, 
and we have only kept the most relevant terms. The 
equation of state is 



We first consider the loss of stability of the symmetric 
(paraelectric, paraelastic) phase at the phase transition. 
According to standard procedure, it corresponds to the 
first appearance of a non-trivial solution to the linearized 
equations of state. We shall look for the non-trivial so- 
lution for the cc-component of the displacement vector 
u x = u(y, z). One can swap the x and y axes to consider 
the u y strains instead. This will only be the considera- 
tion of the inhomogeneous part of the strain. The homo- 
geneous strain of the whole sample defines the change of 
its volume and shape and is described by six independent 
components of the (homogeneous) strain tensor. In our 
case, there is no loss of stability with respect to homoge- 
neous deformation since it would cost an infinite elastic 
energy to produce such a strain with an infinitely thick 
substrate. 

To find the inhomogeneous part of the strain u x at the 
phase transition one should satisfy the equations of local 
equilibrium, doi k jdx k = 0, which in the present case 
read 



da xy da xz 



= 0. 



dy dz 
We shall use the Fourier representation 



u(y,z) 



Uk (z) exp (iky) dk 



(3) 



(4) 



and find the first appearance of the non-trivial solution 
for u for a given wave vector k. We then determine the 
k where the instability sets in first, and this will be the 
point of the stability loss of the symmetric phase. 

We obtain the following equations for the strain with 
the use of Eqs.(|),(|) 



d 2 u k A k , 



dz 2 



—k 2 u k = 0, 0<z<l; 



a xy = 2{A-DV 2 )u 



(2) 



d 2 u k 



k u k — 0, — oo < z < 0, 



(5) 
(6) 
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FIG. 1. Schematic of the domain structure (top panel) with 
a period 2a in a thin ferroelastic film of thickness I attached 
to a substrate (bottom panel). The relative displacement of 
the domain walls by S would change the homogeneous strain 
in the bulk of the film. Since the appearance of the homo- 
geneous deformation in the film is very costly, the change in 
the domain width 5 in the field will be suppressed, and the 
dielectric response would not be enhanced in comparison with 
a monodomain free-standing film. 

where Ak — A + Dk 2 . At the free surface (z = I) the 
boundary condition reads a xz (I) = 0, which is equiva- 
lent to duk{z)/dz = 0. In addition, the displacement 
Uk (z) and the stress o~ xz (z) should be continuous at the 
interface z = 0, and the stress should vanish at z — ► — oo. 

Let us first consider the case of Af. < 0, which would 
lead to a loss of stability of the paraphase. The solution 

ofEqS. §), (|) is 

Uk(z) — F cos r]k(z — l) 7 < z < I; (7) 
Uk{z) = Gexp \k\z, — oo < z < 0, (8) 

where rj 2 = —Ak/fi. The boundary conditions give us the 
condition of the existence of the non-trivial solutions 



cot rjkl = r). 



(9) 



This equation does have a solution for Ak < 0, while 
there is no solution at Ak > 0, hence the loss of stability 
takes place for Ak < 0. For the region of interest (77 <C 1) 
the approximate solution is 77 ~ ir/2kl or 



A I -Dk 2 = 



Ak 2 l 2 

The minimum value | A | c corresponds to 
m [ 2> d i/2 iV2 



(10) 



(11) 



1 /2 

where usually (not on very "soft" substrates) {D/fx) ' ~ 
d a t, with d a t the interatomic distance (see, e.g. H). This 
means that the loss of stability of the symmetric phase 
takes place at 



A 



itD 1 ' 2 ^' 2 



dat 



or at 



T = T cl =T c - 



7rDVV/a 



a/ 



(12) 



(13) 



The temperature of the phase transition is lowered by 
about T at j, where T at is the characteristic "atomic" tem- 
perature. Since normally T at ~ (10 2 — 10 3 ) T c , one may 
expect a complete suppression of the second order tran- 
sition for films with thicknesses of hundreds of atomic 
layers depending on the particular materials parameters. 
The inhomogeneous structure that forms when the sys- 
tem looses stability is a domain structure close to the 
phase transition, where the widths of the domain walls 
and the domains themselves are comparable (see, e.g. 

||). One can readily see that this is valid in the present 

r /2 

case since the domain wall width W = (D/2 \ A |) . 

We consider next the domain structure not very close 
to the phase transition in a state with the spontaneous 
strain u®,. There the domain wall width is much smaller 

y 

than the width of the domains and one can use the lin- 
earized equation of state, obtained by expanding the free 
energy ([!]) about the spontaneous deformation, 



T xy = 2M(u xy -u xv ), 
i° xy = uq — + (~A/B) 1/2 . 



(14) 
(15) 



where M = —2A is -C A* when the system is close to 
the transition ("soft" modulus), with the gradient term 
being the origin of the domain wall energy. 

Initially we shall assume that all the domains have the 
same width, which we will find by minimizing the sum of 
the elastic energy and the (surface) energy of the domain 
walls. We shall follow the standard procedure (see e.g. 
Q). Firstly, we consider the film without a contact to 
the substrate where we create a distribution of the spon- 
taneous strain (a domain structure) in such a way that 
there will be no stresses in the film. The total length of 
the contact area of the film remains exactly the same as 
if it were in the symmetric phase, Fig. 1. Then we join 
the film and the substrate, maintaining the continuity 
of the displacements at the interface. After the contact 
is made, the inhomogeneous stresses appear both in the 
film and in the substrate, but no uniform stresses. 

We consider a stripe-like domain structure with the 
spontaneous strain 

u xy (y,z) — uo, < z < I, (2n — 1) a < y < 2na; (16) 
u xy (yj z ) = ~ u o, < z < I, 2na < y < (2n + 1) a, 

with the period 2a. There would be no stresses in the 
structure if Uq = —A/B. We have to find the displace- 
ments appearing after the film is attached to the sub- 
strate, u x (y, z) = u (y, z). For the film (0 < z < I) Eq.(|^) 
takes the form 



2 



ah 



\f d2u ^2M d2u ° y 
dy 2 ^ dz 2 dy 2 



(17) 



Since the domain pattern is periodic, the elastic displace- 
ments may be represented as a Fourier series 

u (y> z ) = XI Uk ( z ) ex P ( ik y) ' fc = — , (is) 

k 

where n = ±1,±2, . . . After solving the resulting ordi- 
nary differential equations with the above conditions one 
finds the elastic energy by e.g. using the formula 



Fpj 



(19) 



with the result for the elastic (stray) energy per unit area 
of the film: 



Fstray _ 16/lU§a 



A 



-5 — v y 

7T 3 2.7 + 1 



1 



^ (2j + l) 3 cotb^fcZ + T?' 



(20) 



where 77 = yj M//J,. We shall first consider temperatures 
not very far from the transition, where 77 -C 1. Since 
cotha: > 1 for any argument, one can then omit the sec- 
ond term in the dominator. We shall suppose that the 
equilibrium period of the domain structure satisfies the 
condition ijkl 3> 1 and check later that this condition is 
fulfilled. We can then put coth?7fcZ ~ 1 for all terms in 
( pp| ) and find the energy of the domain structure 



F 
A 



14C(3) 2 
5— Liv-oW, 



(21) 



where the first term is the energy of the domain walls, 
with 7 the domain wall surface energy. Therefore, the 
equilibrium period of the domain structure is 



7/ 



14C(3) fiu 2 n 



1/2 



/ 4tt 3 D l /H 
V21C(3) H 1 ' 2 



1/2 



(22) 



since 7 = 8V2D 1 / 2 | A | 3 / 2 jZB |J. We must check now 
if the above assumption Trr/l/a 3> 1 is fulfilled together 
with 77 <C 1. Together those constraints read 



1/2 



.4 



< 1 



(23) 



Since usually (Df/i) 



1/2 



d a t I, the condition is satis- 



fied not very close to the transition. 

The period of the domain structure becomes exponen- 
tially large when the modulus /x is anomalously small, 

(D / n) 1 / 2 ^> d a t, and the condition (|25| ) is violated. Then 

1/2 

for very thin films / <C (D/fj,) one has to reconsider 
the calculations of the sum in Eq. (|2^) . The exact result 
for M = fi (77 = 1) is 



AF 



stray 



A 



; C(3)-Lz 3 ( e - b ) + -Lz 3 (e 



-2b\ 



(24) 



where b = 2irl/a, Li n (z) = Efeli ^A" §, cf - Refs - 
P,^0[. The period of the domain structure for a very 
thin film on a "soft substrate" is exponentially large (cf. 
Refs. §10) ), 

27; = 1 ^exp( TrT7¥7 ). (25) 



7rZ / 7T7 

a= ^ exp U^ 



Next, we shall show that there is no instability with 
respect to the structure with opposite domains of un- 
equal widths, as suggested by Roytburd 0. To this end 
we shall estimate the change of the energy of the system 
where the domains have different widths, a(l + 5) and 
a(l — 6). As we have discussed above, the homogeneous 
part of the strain is zero. This is possible with domains of 
unequal width only if the homogeneous strains in the do- 
mains, u\ = uq + Aui and u% = — uq + A.U2, are different. 
The condition of zero homogeneous strain reads 



Ui(l + <S)+U2(l-<S)=0. 



(26) 



We must also minimize the elastic energy of this struc- 
ture. The energy density of the homogeneous deforma- 
tion can be found from the energy density in the fer- 
rophase, 2M(u xy — u xy ) 2 + 2^iu 2 xz [cf. Eq.(Q)] which 
gives us the elastic energy of homogeneous stresses (cf. 
Ref. §) 

Fh ' ' 1 — <^ 2 1(1 -S)2M{Au 2 ) 2 . (27) 



A 



i(l + *)2M (A Ul ) 2 



Minimizing Eq.(p7|) with respect to Aiti (which is equiv- 
alent to the condition that the force acting on the domain 
wall is zero), taking into account the constraint (p6|), we 
obtain Aui = Ai*2 = —uq5 and 



F h /A 



2Mu 2 lS 2 . 



(28) 



One has to add the energy of inhomogeneous stresses, 
which depends on 5 



AF, 



stray 



A 



8fiularj ^ (-1)™ (1 - cos nnS) 

^3 1^—7 — -. I 29 ) 



n=l TV 



(coth 



1 



where AF stray (8) = F stray (S) - F stray (0) , to obtain the 
total change in energy. 

We shall now analyze two limiting cases. 
Rigid substrate not very far from T c (77 = y/M/jj, <c 
1, irql/a 3> 1)-— The expression in parentheses in the 
denominator in Eq. ( p9[ ) can be replaced by unity, and we 
obtain 



AF 



stray 



8/i7jQa77 



A 



7T° 

4 In 2 

7T 



-C(3)+ReLi 3 (-e i7r5 ) 



fiUQarjS 2 



(30) 



3 



In spite of the negative sign of this stray contribution, 
the total energy increases with 5, 

— = 2 W u 2 a( V )8 2 > 0, (31) 

J\ a 7r 

and, consequently, no instability sets in. 

Soft substrate [{D/fi) 1/2 > d at , r\ = 1].- There the 
contribution (|29|) can be found exactly with the result B] 



AF s tray (S) 

A 



—fiu a 

7T 



r 3 
+Li 3 



Re Lia 



„i7r5— 6> 



(32) 



In this case the domains may be wide, l/a <C 1, Eq.(25) 
and we obtain 



AF 



;2 

) I 



— = 7T^U 

.A a 



2 <S 2 . 



(33) 



This energy increases with S and, therefore, the domain 
structure remains stable on a "soft" substrate too. 

Finally, we shall evaluate the dielectric response of the 
above domain structure. We assume a linear coupling 
between u xy and P z given by the thermodynamic poten- 
tial 

F = J dV{2A 1 u% + \a 2 P 2 + Bu% + gu xy P) (34) 

with A2 > 0, P = P z . Minimizing ([34j) with respect to 
P we obtain 



AoP 



0. 



(35) 



hence the spontaneous strain ±uq leads to the sponta- 
neous polarization in the domains Pq = ^fguo/A2. The 
previous free energy functional ([!]) can be obtained from 
(|J) with the use of (f§), and A = A X - g 2 /4A 2 . 

Let us show now that the net polarization does not 
change when the domain walls have moved, so that now 
(5^0. In a free standing film the change would be 2P S, 
but in a film on a substrate the polarization itself changes 
since it follows the changes of the mean strains in every 
domain. The homogeneous polarization Ph is then, with 
the use of Eq. (ph, 



(36) 



P h = \{l + &)Pi + \{l-5)P2 



2A, 



[(l + 5)u 1 + (l-5)u 2 } = 0, 



because the average strain in the film (square brackets) 
is exactly zero, Eq. (p6|). In other words in a thin film 
on a substrate the shift of the domain walls does not 
lead to any change in the average polarization. They 
appear to be decoupled, while they are strongly coupled 
in a free standing film. Surprisingly, under an external 
field the domain walls do not shift in our case and the 



field produces Ph = E/A2: one has to put the field E 
in Eq.(Q) and set u xy = because the homogeneous 
deformation remains zero. 

The result ([36]) emphasizes the qualitative difference in 
the dielectric response of a ferroelectric-ferroelastic do- 
main structure in a free standing film and in an epitaxial 
film: (i) there is a finite stiffness with respect to the 
domain wall shifts and (ii) these shifts change the polar- 
ization within the domains for the film on a substrate. 
In our specific case these changes completely compen- 
sate the polarization gain that would have taken place in 
a free standing film. In general, the compensation is not 
perfect but the effect should be taken into account. For 
example, these two aspects have been overlooked by Erbil 
et al. |ll[ who have only accounted for the energy change 
when the variation of the relative domain fraction. This 
has been ascribed to a pinning potential, while the other, 
much more important, terms were not considered. Their 
conclusion about enhancement of the permittivity due to 
the domain structure is incorrect. Our results, addition- 
ally, do not show any signs of instability of the domain 
structure with the film thickness, as has been speculated 
first by Roytburd [Q analytically and then by Pertsev 
et al. |Q numerically for perovskite ferroelectrics. Since 
their case is somewhat different (the spontaneous strain 
Uxx — Uzz), it might be interesting to apply their numerics 
to the present simpler case to check the result. 

In conclusion, we have demonstrated analytically how 
the domain structure sets in at the phase transition and 
that it remains stable with respect to spontaneous break- 
ing of symmetry between the opposite domains irrespec- 
tive of whether the substrate is rigid or "soft" . The di- 
electric response of the epitaxial films is qualitatively dif- 
ferent from that of free standing films. In the present 
case the motion of the domain walls has zero effect on 
dielectric response in ferroelastic-ferroelectric thin films. 
Generally, the response of an epitaxial thin film is sup- 
pressed compared to a free film. 
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